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Abstract 

Using a combination of the operator product and heavy quark 
expansions, we resum the leading nonperturbative contributions to 
the inclusive photon spectrum in B — * X s 7 decays. The shape of 
the spectrum is determined by a universal structure function, which 
describes the distribution of the light-cone momentum of the 6-quark 
inside the £>-meson. The moments of this function are proportional 
to forward matrix elements of higher-dimension operators. As a by- 
product, we obtain the bound Ai < for one of the parameters of 
the heavy quark effective theory. The integral over the B — > X s 7 
structure function is related to the shape function that governs the 
fall-off of the lepton spectrum close to the endpoint in B — » X u I v 
decays. A measurement of the photon spectrum in rare -B-decays can 
therefore help to obtain a model-independent determination of V u b- 
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1 Introduction 



The rare decays of B-mesons are of great importance, since they are sensitive 
probes of new physics beyond the standard model (see, e.g., Refs. Jj], |J). Such 
processes are induced by the exchange of heavy particles, which manifests 
itself at low energies in the appearance of local operators multiplied by small 
coefficient functions. These coefficients depend on the masses and quantum 
numbers of the heavy particles. In the standard model, rare decays of the 
type b — > 57 are induced by penguin diagrams with virtual top or charm 
quark exchange. Recently, the first such decay mode, B —>■ K*j, has been 
observed by the CLEO collaboration. The reported branching ratio is (4.5 ± 
1.9 ± 0.9)% 0. The interpretation of this result is difficult, however, due to 
the lack of reliable theoretical methods to calculate the low-energy hadronic 
matrix element of the penguin operator between meson states. Existing 
estimates of this matrix element rely on quark models @, || or QCD sum 
rules m 0, §] and are thus model-dependent. 

For several reasons, one expects that the theoretical analysis is more 
reliable for inclusive B — > X s 7 decays, where one sums over all possible final 
states containing a strange particle. Assuming quark-hadron duality, the 
inclusive decay rate was traditionally calculated using the free quark decay 
model and including short-distance corrections from virtual and real gluons 
[|, 0, [LI], II, m, m, U, U, |l7|. Recently, however, it has been observed that 



inclusive decays of hadrons containing a heavy quark Q allow for a systematic 
expansion in powers of A/rriQ, where A is a characteristic low energy scale 
of the strong interactions [TjJ |2D|, [2^, |23|. The parton model emerges 
as the leading term in this QCD-based expansion, and the nonperturbative 
corrections to it are suppressed by a factor K 2 /m 2 Q. The fact that there are 
no first-order power corrections relies on a particular definition of mg, which 
is provided in a natural way by requiring that there be no residual mass 



term for the heavy quark in the heavy quark effective theory [24, B^]. This 
definition is unique and can be regarded as a nonperturbative generalization 
of the concept of a pole mass. 

The availability of a systematic expansion of the inclusive B — ► X s 7 
decay rate raises the hope for a better understanding of rare decays, which 
is necessary to increase the sensitivity to new physics. For practical reasons, 
however, it is not sufficient to have a reliable calculation of the total decay 
rate. In fact, the distribution of the photon energy will be affected by various 
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experimental cuts, and it is thus the spectrum dT/dE^ that needs to be 
calculated. In the free quark decay model, the photon in b — > 57 decays is 
monochromatic. Corrections to this simple picture arise from two sources: 
Real gluon emission produces three-body final states, leading to a continuous 
energy spectrum. These effects have been calculated in perturbation theory 
nj; they will not be discussed here. In addition, bound-state corrections in 



the initial state, in particular the "Fermi motion" of the 6-quark, lead to a 
dispersion of the spectrum. These nonperturbative effects can Doppler-shift 
the spectrum above the parton model endpoint. So far, such effects have 
been estimated |26| using the phenomenological model of Altarelli et al. . 

In this paper, we present a more rigorous treatment of bound-state cor- 
rections to the free quark decay model. In particular, we show that QCD 
provides a natural framework to account for the "Fermi motion". Extend- 
ing our previous analysis of the endpoint region of the lepton spectrum in 
B — > X u tv decays [ffiH , we resum the leading nonperturbative contributions 
to the photon spectrum in B — > X s 7 decays to all orders in the 1 / m& expan- 
sion. We show that the spectrum is determined by a fundamental structure 
function, which describes the light-cone residual momentum distribution of 
the heavy quark inside the S-meson. Quite remarkably, an integral over this 
structure function is related to the shape function that governs the endpoint 
region of the lepton spectrum in B — > X u t v decays. In Sect. ^|, we generalize 
the results of Refs. [19|, ^2] and construct the operator product expansion for 
the photon spectrum in B — > X s 7 decays, including the leading nonpertur- 
bative corrections. In Sect. |3|, we perform a resummation of the most singular 
terms in the photon spectrum to all orders in I/ttv The moments of the 
spectrum are related to forward matrix elements of local, higher-dimension 
operators in the heavy quark effective theory. We show that the characteris- 
tic width of the spectrum is determined by the kinetic energy of the 6-quark 
inside the I?-meson. From this relation, we derive the bound Ai < for 
one of the fundamental parameters of the effective theory. In Sect. [|, we 
illustrate our results using a toy model, which is a simplified version of the 
approach of Ref. P7| . In Sect. we derive the relation between the photon 
spectrum and the 6-quark structure function of the S-meson, and we relate 
the Fourier transform of the spectrum to the forward matrix element of a 
gauge-invariant, bilocal operator. Section [|is devoted to a discussion of the 
connection between rare -B-decays and inclusive semileptonic decays. We 
show that the shape function ||28|, which describes the fall-off of the lepton 
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spectrum close to the endpoint in B — > X u £u decays, is given by an inte- 
gral over the photon spectrum in rare B — > X s ^ decays, up to corrections 
of order l/m b . This connection may help to obtain a model-independent 
determination of V ub . In Sect. [7|, we summarize our results and give some 
conclusions. 



2 Operator Product Expansion 

In this section, we discuss the application of the operator product expansion 
to the inclusive rare decays B — > X s 7. We derive expressions for the pho- 
ton spectrum and the total decay rate, to order l/m| in the heavy quark 
expansion and to leading logarithmic order in a s (m(,). In the limit where the 
mass of the strange quark is neglected, the total decay rate and the average 
photon energy have been calculated in Refs. p^j . We will generalize the 
results presented there. 

In leading logarithmic approximation, the rare decays of interest are me- 
diated by an effective Hamiltonian containing a local penguin operator, 
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H cS = -k V tb V t l c 7 (m b ) — s <j*» (m b P R + m s P L )b F„ v , 



where Pl = |(1— 75) and Pr = §(1+75) are left- and right-handed projection 
operators, and is the electromagnetic field strength tensor. The Wilson 
coefficient c 7 {rrib) describes the evolution from high-energy scales \x ~ m t or 
m w to low-energy scales /1 ~ m b [g |TD|, [11], [TJ, [TJ, [T|, |IJ, [lj| 0. It is 



sensitive to the mass of the top quark and, more generally, to any kind of 
new physics beyond the standard model. 

In terms of the effective Hamiltonian, the inclusive differential decay rate 
is given by 

dr= ^9 d ?9 F E (2vr) 4 5V^-P7-^)|(X s (^) 7 (p 7 )|^e ff |^)) 

(2) 

where we sum over the two transverse polarization states of the photon. We 
use a mass-independent normalization of states such that 

(B(v)\B(v))=v (2n) 3 6 3 (0), (3) 



3 



where v is the four-velocity of the 5-meson. The decay rate can be written 
in terms of the imaginary part of a correlator of two local currents, which 
contains all dependence on hadronic dynamics. We define 



T(v,p 7 ) = -i I dx e i{ - mbV - p ^ x (B(v)\ T{J»(x),Jl(0)}\B(v)) , (4) 

where 

J"(z) = b v (x) [Y,^] (m b P L + m s P R ) s(x) ; b v (x) = e im ^ x b(x) . (5) 



This leads to 



d% G 



dF = T^hr ^TT I V »K\ 2 \c 7 (m b )\ 2 lmT(v,p,) . (6) 



(2tt) 3 2E 7 8 
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Note that the Fourier components of the rescaled heavy quark field b v (x) in 
(HD contain the "residual" momentum k = p b — m b v. In contrast with the 
full heavy quark momentum p&, the residual momentum is of order Aqcd- It 
is then appropriate to construct an expansion in powers of k/rrib- 
To this end, it is convenient to introduce dimensionless variables 

P-y r- m s 2i? 7 

P= — , m = VP= — , y = 2v-p= , (7) 

777,6 "Mb TTlb 

where denotes the photon energy in the rest frame of the S-meson. Since 
v 2 = 1 and p 2 = 0, the function T(v,p^) only depends on the kinematic 
variable y. The correlator T(y) is analytic in the complex y-plane, with 
discontinuities on the real axis. The physical region corresponding to the 
decay B — > X s 7 is 

< s <!2»(i-!*Wo9, (8) 

where we have assumed m b ~ 4.7 GeV for the purpose of illustration. In 
addition, there is a cut starting at y ~ 4.16 corresponding to the process 
7 + B — > X s bb, where X sbb contains two 6-quarks and an s-quark. This 
unphysical cut is separated from the physical one by a large energy gap 
AE^ = ms (1 + mK*/rriB) 2 — 7.2 GeV. Because of this analytic structure, 
phase-space integrals of T(y) with smooth weight functions can be deformed 
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from the physical region into contour integrals far away from the physical sin- 
gularities. It is then possible to construct an operator product expansion of 
the correlator in terms of local operators Cj, which contain the 6-quark fields 
and have dimension d > 3 |T8|, [19|, |20|, |21], |22|, |23| . These operators are mul- 



m, 0, 123, 

tiplied by short-distance coefficient functions Ci(y), which can be computed 
in perturbation theory using free quark states. The leading contributions 
in a s come from a tree diagram with an intermediate s-quark carrying the 
momentum (pb — p 7 ) = (m&i; + k — p 7 ) . The operator product expansion 
is obtained by replacing the residual momentum A; by a covariant derivative 
iD. This gives the propagator in the background field of the light degrees 
of freedom in the decaying 5-meson. Next, one expands the propagator in 
powers of iD/rrib. This gives 



m b($ — § — rh + it) + ip 



f — § + rh 
nib A 



E 



n=0 L 



ip 



{i> —ff+rh) 
nib A 



where 



A 



p) 2 — rh 2 + ie = 1 — y — p + ie . 



(9) 



(10) 



Once the Wilson coefficients have been determined, one has to evaluate 
the forward matrix elements of the local operators Oi between 5-meson 
states. The leading operators have dimension three and can be related to 
matrix elements of vector and axial vector currents. The vector current 
matrix element is normalized by current conservation, 



(B{v)\l vl »b v \B{v)) = rf 



(11) 



whereas the matrix element of the axial vector current vanishes by parity 
invariance. These leading-order contributions reproduce the parton model. 
The nonperturbative corrections to it are described by the matrix elements 
of higher-dimension operators. They can be evaluated using the powerful 
formalism of the heavy quark effective theory The heavy quark field b v 
is split into "large" and "small" two-component spinors h v = |(1 + f) b v and 
H v = |(1 — f) b v , and the field H v is integrated out to obtain an effective 
Lagrangian pOl 31L |32|. Any operator in the full theory has an expansion 



in terms of operators in the effective theory, which only contain the field 
h v . Using techniques |22], [HJ that are standard by now, we obtain the 
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relations 



(B(v) | l v r iD» b v \B(v)) = Tr {r (y + ^ - 5<AP + )} , 

6 ' 

+ ^Tr{TP + ia^P + }, (12) 

where P + = |(1 + ^), and T denotes an arbitrary combination of Dirac 
matrices. The low-energy parameters Ai and A2 are related to the kinetic 
energy K b of the heavy quark inside the 5-meson, and to the hyperfine mass 
splitting between B- and f?*-mesons. They are defined as 

K b =-^-, m|.-m| = 4A 2 . (13) 

Using these results, we find 

r(y) = _4 m g( 1 + P )^| 1 _^/V + 7y +5 



A I 6m 2 b \A 2 A 

A 2 f5y-8 



2m\ V A 



5)+ OK" 3 )}, (14) 



with A as given in fllPf ). In the limit p = 0, this agrees with Ref. p2| . For 
the inclusive photon spectrum, we obtain from (^|) 

f = I V M 2 M™>)l 2 (1 + P )(l - p) 3 % s(y, p) , (15) 

where 

Ai — 9k\ 2 3 + 5p 

77 & = 1 + \ 2 + OK" 3 ) ; « = — — ^ . (16) 



The spectral function s(y,p) is given by 



p) = 5(1 - y - p) - (1 - p) 1 1 2 2 - y - p) 



Ai + 3/tA2 s-,, 
( 

(l-rf 2 Af (1 _^^0 K 3). (17) 
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The interpretation of the singular structure of this function is the main sub- 
ject of this paper. Recall that the operator product expansion, which led 
to QT5J) , was only justified when the spectrum is integrated with a smooth 
weight function. Hence, one should understand the singular expression (|I~7j) 
in the sense of distributions. Integrated quantities such as the total decay 
rate and the average photon energy obey a well-defined l/m& expansion. We 
find 



i 



32tt 



(18) 



In the limit p = 0, we confirm the results of Refs. [[HJ |2"2]| . 

In order to obtain an estimate of the magnitude of the nonperturbative 
corrections, we use the quark masses m& = 4.7 GeV and m s = 0.2 GeV, 
corresponding to p ~ 2 x 10 -3 . From the observed value of the B*-B mass 
splitting, one obtains A 2 — 0.12 GeV 2 . The parameter Ai is not directly 
related to an observable. The field-theory analogue of the virial theorem 
relates the kinetic energy of a heavy quark inside a hadron (and thus Ai) 



to a matrix element of the gluon field strength tensor [5^]. This theorem 
makes explicit an "intrinsic smallness" of Ai, which was not taken into ac- 
count in existing QCD sum rule calculations of this parameter |RSL |3"7| ]. 
As a consequence, we expect that (— Ai) is smaller than predicted in these 
analyses. Here we shall use the range — Ai = 0.1-0.3 GeV 2 . We then obtain 
rjb ~ 0.97, corresponding to a 3% decrease of the parton model decay rate. 
The correction to the average photon energy is below 1%. 



3 Resummation of the Most Singular Terms 

The numerical analysis presented above shows that the nonperturbative cor- 
rections to integrated quantities such as the total decay rate are very small. 
Were it just for these corrections, one could say that the most important 
result of the formalism presented in the previous section would be to give 
a theoretical justification for the parton model. However, in this section we 
will show that much more interesting information can be extracted from this 
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analysis. This information is encoded in the coefficients of the singular terms 
in the spectral function s(y,p). 

The singularities at y = 1 — p in flTTD arise from the fact that the operator 
product expansion becomes singular when the s-quark propagator is almost 
on-shell f22|. Although it is legitimate to evaluate phase-space averages of 



smooth functions using the singular theoretical expression, one cannot trust 
the shape of the spectrum as given in fll7D . The true spectrum will be differ- 
ent. A similar situation is encountered in inclusive semileptonic B — > X u iv 
decays, where the lepton spectrum obtained from the operator product ex- 
pansion becomes singular when the lepton energy approaches the parton 
model endpoint |19|, |21]. In Ref. [p8|l , we have suggested that the devi- 



ation of the true lepton spectrum from the parton model prediction could 
be described by a shape function S(y), which has a support only in a small 
region close to the endpoint. The singularities in the theoretical spectrum 
can be identified with the first few terms in a moment expansion of the shape 
function. 

We can adopt a similar point of view in the case of B — > X s 7 decays. Let 
us identify the function s(y,p) in (|T5| ) with the physical photon spectrum 



subject to the normalization condition 

00 

Jdys(y,p) = l, (19) 


so that the total rate is given by the first equation in (|T8|). We will show below 
that the characteristic width a y of the spectrum is proportional to l/m&. 
When integrated with a smooth function that is slowly varying on scales 
of order l/m b , the spectral function s(y,p) can be replaced by a singular 
expansion [28 



s(y, P ) = f^^\l-y-p). (20) 

n=0 U - 

The moments M n (p) are defined as 

00 

M n (p) = Jdy(y-l + p) n s(y,p). (21) 


By definition, M (p) = 1. We can now identify the singular expression in 
(|T7D with the first few terms in this expansion and relate the moments M n (p) 
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to nonperturbative parameters. Neglecting terms of order l/m|, we obtain 



2m{ 



M 2 (p) = -(l-p) 2 A^. (22) 



From dimensional analysis, it follows that the moments obey an expansion 
of the form 

M„W = ^ + ^# + ..., (23) 



m h 



with coefficients a n (p) and b n (p) that are independent of rrib (up to logarithms 
arising from radiative corrections). Hence, the QCD prediction that Mi(p) 
is of order 1 /ml indicates a non-trivial cancellation: The shift in the average 
value of y due to bound-state effects is of order 1 /ml, corresponding to a 
shift of order AQ CD /m;, in the average photon energy. Naively, one would 
expect this shift to be of order Aqcd- 

The second moment is a measure of the width of the photon spectrum. 
As stated above, we find that a y oc X/m^. In Ref. |[28|| , we showed that the 
quantity a y evaluated for p = describes the width of the endpoint region 
of the lepton spectrum in B — > X u i v decays. The connection between these 
two cases will be discussed in more detail in Sect. |[ Note that the width 
is determined by the parameter (— Ai), which is proportional to the kinetic 
energy of the 6-quark inside the S-meson [see (|X3|)] . Since, by definition, the 
second moment is positive, we obtain the upper bound Ai < 0. Although this 
result is certainly not surprising, it is not trivial, since operator renormaliza- 
tions could spoil the positive definiteness of the kinetic energy operator that 
defines (-Ai) ||. 

For an understanding of the properties of the spectrum close to the end- 
point, it is not sufficient to truncate the l/m^ expansion at order 1/ml [E9]. 
What is relevant are the rescaled moments 



7d£ 7 [2£ 7 - m h (1 - p)] n s(E 7 , p) = a n (p) + + . . . , (24) 

J J m\) 



where s(E ry ,p) = (2/m b ) s(y, p), such that / °°d£/ 7 s(E^, p) = 1. These mo- 
ments are all equally important in the limit m& — > oo. Hence, it is necessary 
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to resum the operator product expansion. However, with the exception of the 
first moment, for which the coefficient a\{p) vanishes, one may argue that it is 
a good approximation to keep the leading coefficient a n (p) for each moment. 
The corrections involving b n (p) change the moments by small amounts. 

The aim is thus to construct a partial resummation of the operator prod- 
uct expansion, in which one keeps the leading term in each moment but 
neglects l/m& corrections. A crucial observation is that, at any order in the 
1/mb expansion, the coefficient a n (p) receives contributions only from the 
most singular terms in the theoretical expression for s(y,p). For instance, 
at order I /ml the coefficient of the ^"-function in (pi) determines 0,2 (p)', the 
coefficient of the (^'-function, however, determines foi(p). A resummation of 
the most singular terms can be constructed by using the following alternative 
way of writing the s-quark propagator in (0): 



f — ff+ rh 



m b {j) — § — m + ie) + ip m b A + 2(v — p) ■ iD 

1 



00 

n=0 L 



(j) — § — m) ip 



m^A + 2(v — p) ■ iD 



(25) 



Note that all terms but the first one are multiplied by a factor A = {jj — 0+ 
m){f — fj — rh). Since A vanishes at the endpoint, it follows that the first 
term is more singular than the other ones. Hence, the leading singularities 
can be resummed by using the replacement 



fj '+ m 



vribiji — fj ' — rh + ie) + ip m^A + 2(v — p) ■ iD 



(26) 



for the s-quark propagator. The imaginary part of this expression is given 
by a 5-function, and it is straightforward to find that 



s(y,p) 



1 — y — p H (v — p) • iD 



+ less singular terms, (27) 



where we define the expectation value of an operator O as 

(B(v)\h v Oh v \B(v)} 



O) 



(B(v)\h v h v \B{v)) 



(28) 



Here, h v is the velocity-dependent heavy quark field in the heavy quark effec- 
tive theory |52"fl , and the states are the eigenstates of the corresponding 
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effective Lagrangian.[] Equation ( p7[ ) is a formal definition of the spectral 
function, which is valid to all orders in the l/m?, expansion. The "less sin- 
gular terms" omitted here do not contribute to the leading-order coefficients 
in the expansion of the moments. Expanding our result in powers of l/m&, 
and comparing with (0), we find 



a n {p) = ( [2{v-p) -iD] n ) 



y=l~p 



(29) 



In order to extract further information from this relation, it is necessary 
to investigate the structure of forward matrix elements of higher-dimension 
operators in the heavy quark effective theory. Using the equation of motion, 
iv • D h r 



0, one can show that |28 

(iD») = 0, 
{iD»iD y ) = A 2 (/«" 

( iD* iD u iD a ) = A 3 - g^ a ) v l 



(iD^ 1 ... iD^) 



A n v^ 



+ terms with 



l-h 



(30) 



where A2 = — Ai/3 |Jv|]. For n > 4, the matrix elements can no longer 
be parameterized by a single parameter A n . However, terms involving the 
metric tensor give only small contributions, of order p, to the coefficients 
a n (p) in (|29|). Hence, we obtain a (p) = 1, ai(p) = 0, as well as 



2 -^1 



a 2 (p) = -(1 - P) 
a 3 (p) = (l-p) 2 (l+p)A 3 , 
a n {p) = A n + 0(p); n>4. 



(31) 



Let us summarize the main results of this section: Apart from radiative 
corrections, the inclusive photon spectrum in B — > X s ~f decays can be de- 
scribed by a spectral function s(y,p), which is a genuinely nonperturbative 
form factor that accounts for bound-state effects in the decaying meson. The 
moments of this function obey a well-defined 1 /m& expansion. The leading 

1 It is sufficient to evaluate the matrix elements in the effective theory, since we are not 
interested in the higher-order corrections in (|25|). 
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terms in this expansion are related to forward matrix elements of higher- 
dimension operators in the heavy quark effective theory. In the limit p = 0, 
these matrix elements are described by a set of fundamental parameters A n . 
Since the even moments of the spectral function are positive, it follows that 
A 2n > 0. In particular, this gives the upper bound Ai < 0. 



4 A Toy Model 

Before we discuss in more detail the physics of the results derived in the 
previous section, we find it instructive to illustrate them in the framework of 
a toy model, which is a simplified version of the phenomenological approach 
of Altarelli et al. (ACM) [0. In the ACM model, the validity of the parton 
model is assumed, and bound-state effects are incorporated by assigning a 
momentum distribution </>(|pb|) to the heavy quark.0 It is then appropriate to 
replace the covariant derivative by the spatial components of the heavy quark 
momentum p&. The gluon field and the time-component of the covariant 
derivative are neglected. Accordingly, in the rest frame of the 5-meson, one 
makes the replacement 

2{v-p)-iD (32) 

y=l-p 

where p» = pb ■ p 7 /|p 7 | denotes the component of the 6-quark momentum in 
the photon direction. The matrix elements in (B^) and (|29| ) are replaced by 
integrals over the momentum distribution of the heavy quark. In the ACM 
model, one assumes a Gaussian momentum distribution, 



0(N) = 7^i ex H-7rJ' (33) 

where p-p is the Fermi momentum. It is straightforward to calculate the 
hadronic matrix elements a n (p) in this toy model. We find 

w(^-l)!! n 
ill — p) 7s — Vv ! n even, 

^(p) = (l-p) n (pV = { ^ (34) 

; n odd. 



2 In addition, the heavy quark mass is treated as a momentum-dependent parameter 
rrib(\pb\). For simplicity, we shall not consider this aspect of the model. 
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Comparing ( |34] ) with ([H]) for n = 2, one obtains the relation — Ai = |pp 
between the low-energy parameter Ai and the Fermi momentum [2£|. Finally, 



it is possible to calculate the leading term in the spectral function s(y,p) in 
The result is again a Gaussian distribution: 

s, oy (y,p) = ■jtL- exp { - {y -^ p) } . (35) 



The width a y has been defined in fl22|). Such a Gaussian model for the spectral 
function, supplemented by a radiative tail extending to smaller values of y, 



has been used in the analysis of Ali and Greub ||26|| . 

At this point, we have to stress that our toy model is presented for peda- 
gogical purposes only; we do not claim that it provides a realistic description 
of the spectral function. In fact, this simple model of the "Fermi motion" is 
inconsistent with QCD. Note that the vanishing of the odd moments of the 
distribution function, which implies the symmetry of the spectral function 
around y = 1 — p, is a consequence of rotational invariance. As such, it is 
unavoidable in a model where the time-component of the heavy quark mo- 
mentum is neglected. In QCD, however, there is no reason why any of the 
coefficients a n (except ai) should vanish. Hence, we expect that the physical 
spectral function is asymmetric. In the following section, we shall discuss 
in more detail the correct generalization of the toy model in the context of 
QCD. This will lead us to the concept of a universal light-cone structure 
function, which replaces the distribution function 0(|p;,|). 



5 The Light-Cone Structure Function 

The alert reader will have realized the close analogy of our discussion in 
Sect. [5] with deep inelastic scattering. In this section, we will exploit this 
relationship. From now on, we will neglect the mass of the strange quark 
and set p = 0. We expect this to be an excellent approximation. For instance, 
the coefficients ci2(p) and a^{p) in (|3lD change by less than 1.5% when m s is 
varied between and 0.4 GeV. 
For p = 0, the vector 

rif, = 2(v-p) fl (36) 
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is a null vector on the forward light cone satisfying n 2 = and n ■ v = 1. Let 
us denote the scalar product of a four- vector p with n by n • p = p+. In the 
rest frame of the 5-meson, we are free to choose n M = (1, 0, 0, 1), such that 
P+ = P° + P 3 - Moreover, we can simplify expressions by using the light-cone 
gauge (LCG) n ■ A = 0. From (gD and ([31]), it then follows that 

LCG 

a n (0)=A n = ((iD+) n ) = ((id + T). (37) 



This is the correct generalization of (^)- Note that id + is the operator 
corresponding to the light-cone residual momentum k + of the 6-quark in the 
I?-meson. 

Using this notation, the spectral function s(y) = s(y, 0) in ( p7|) takes the 
form 

s(y) = ~~ V + " — + l ess singular terms 

= Jdk + 5^1-y+^j{f(k + ) + 0(l/m h )}, (38) 

where 

f{k + ) = {d(iD + -k + )} = - 39 

(5(17)1 ft„ /i„ \B(v)) 

is a universal structure function, which determines the probability to find a b- 
quark with light-cone residual momentum k + inside the 5-meson. Since this 
function is defined in terms of a matrix element in the heavy quark effective 
theory, it is independent of the 6-quark mass. The moments of f(k + ) are 
given directly in terms of the hadronic matrix elements A n defined in (|30D, 



A n = J dk+kl f(k + ). (40) 

The corresponding moment expansion reads 

00 (—1 ) n \ A 

f(k + ) = J2 [ -^A n 5^(k + ) = 6(k + )-^5"(k + )--^6"'(k + ) + .... (41) 
n\ 6 6 

Equation ( |38|) shows very nicely the structure of the 1/nib expansion 
proposed in this paper: The photon spectrum is a convolution of the free 
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quark decay spectrum (i.e. a (^-function) with a nonperturbative distribution 
function. The leading term in the l/m& expansion of this function is given by 
the universal structure function f(k + ). This is the correct generalization of 
the toy model discussed in the previous section. The important distinction 
is that the light-cone momentum k + contains the time-component of the 
residual momentum k. Thus, the odd moments of the structure function are 
not forced to vanish by rotational invariance, and f{k + ) is in general not 
symmetric around k + = 0. 

One can consider (|38D as the recipe for a systematic and consistent imple- 
mentation of the leading bound-state corrections, even when one goes beyond 
the leading order in perturbation theory. Needless to say, however, the struc- 
ture function f(k + ) cannot be calculated from first principles. One option is 
to extract this universal function from experimental data. The photon spec- 
trum in rare B — > X s ^ decays is an ideal place for this. Alternatively, one 
may try to obtain a QCD-based prediction for f{k + ) using nonperturbative 
techniques such as lattice gauge theory or QCD sum rules. For this purpose, 
it may be useful to relate the structure function to a forward matrix element 
of a gauge-invariant, bilocal operator. Let us define the Fourier transform of 
f(k+) as 

f(t) = f dk + f(k + ) e- ik + l . (42) 



It then follows that 

z 



(B(v)\h v {0)P exp ( -ifdx^(x)) h v (z) \B(v)) 
i 2 L (43) 

where z = tn is a four- vector on the light cone satisfying z 2 = and z-v = t, 
"P" denotes path ordering, and the integral is along a straight line. In light- 
cone gauge, the phase factor equals unity. The function f(t) describes the 
spatial distribution of the 6-quark inside the 5-meson. 

Finally, we note that the residual momentum structure function f(k + ) 
obeys a simple relation to the usual structure function &b(x), which deter- 
mines the probability to find in the 5-meson a 6-quark with total light-cone 
momentum fraction x. Using that pb = m^v + k and n ■ v — 1, we have 

x,i% = ^, (44) 
[Pb)+ m B 
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and hence 

&b(x) dx — (f(k + ) + 0(l/m&)j dk + ; k + = m B x—m b . (45) 

From the requirement that < x < 1, it follows that the allowed range for 
the light-cone residual momentum is — nib < k + < m B — nib- In the limit 
nib — ► oo, this becomes 

- oo < k+ < A , (46) 

where A denotes the asymptotic value of the mass difference nig — nib and 
can be identified with the effective mass of the light degrees of freedom in the 
.B-meson [^4], [25]. An interesting consequence of ( f46|) is that it determines 
the kinematic endpoint for the photon energy in B — > X s ^ decays. From 
(HD, we find that 

Z/ m ax = 1 + — + 0{nif) = m + 0{m^ 2 ) . (47) 
m b m h 

This is in fact consistent with the physical endpoint given in 

Knowing the moments of f{k+), it is straightforward to calculate the 
moments of the structure function 6b(x). We find 

\ i ra / \ 

/ dx (1 - x) n b B (x) = — E (-!)" (I K A n ~ fc + ^W- 1 ) • (48) 
^ m B fc=o v/cy 

In particular, this leads to the sum rules 

l 

1, 



o 

i 





1 

/ 



J dx b B {x) 
o 

/" A 

/ dx (1 - x) b B (x) = h 0(m„ 2 

J ni B 

dx (1 - xf b B (x) = \ ( A 2 - ^ ) + 0(m B 3 ) . ( 19 ) 



o 

The second relation has been derived previously, in a different context, in 
Refs. [139], E0|. Note that it implies the lower bound A > 0, which in view of 



the recent criticism |38| of the Guralnik-Manohar bound A > 237 MeV []4T 
seems less trivial than one may think. 
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6 Relation to B — > X u £ v Decays 



In this section, we discuss an interesting relation between the leading non- 
perturbative corrections in rare and semileptonic inclusive -B-decays. This 
connection may help to obtain a model-independent determination of the 
element V ub of the Kobayashi-Maskawa matrix. In the limit m u = mi = 0, 
the inclusive lepton spectrum in B — > X u £ v decays can be written as 



dT 

dy 



Gl 



V, 



ub 



96vr 3 



■ml 



F(y)e(l-y)+F(l)S(y) 



y 



2Ee 



(50) 



In this case, the kinematic variable y denotes the rescaled lepton energy, and 
F(y) is a slowly varying function of y. Apart from small nonperturbative 
corrections of order \jm\ [[Tj^ |21| , it is given by parton model kinematics: 
F(y) ~ y 2 (3 — 2y). Close to the endpoint, one can replace F{y) ~ -^(1) — 1, 
up to corrections of order \jm\. The shape function S(y), on the other 



hand, is a rapidly varying, genuinely nonperturbative object |p£;|. It is non- 
zero only in a small region around the endpoint of the spectrum. Using the 
resummation technique developed in Sect. EjL we find that 



9(1 -y) + S(y) 



/ in-D\\ 

B 1 — y -\ ) + less singular terms 

V m b J I 

[ dk + e(l-y + ^-){f(k + ) + 0(l/m b )\. (51) 



Here, n = 2{v — pe/m b ) is again a null vector when y — 1, and hence f{k + ) 
coincides with the universal function defined in ([39]). We observe that, close 
to the endpoint region, the lepton spectrum in semileptonic B —>■ X u I v 
decays can again be written as a convolution of the free quark decay distri- 
bution (i.e. a step function) with the structure function f(k + ). Comparing 
the above relation with (p8|), we obtain 

d r n 

s(y) = ©(I ~ u) + S(y) + l ess singular terms. (52) 

Let us integrate this equation to obtain 



Q(l-y) + S(y) ~ Jdy's(y') 



17 



dy' 0(1 -y>) + S(y') 



dy' (y' - y) s(y') 



(53) 



These relations are exact up to corrections of order l/m^ or a s . 

The main goal of the study of B — > X u £ v decays is to extract the element 
V u b of the Kobayashi-Maskawa matrix. The experimental analysis is very 
complicated, as there is only a tiny window close to the endpoint region where 
the signal is not overshadowed by the large background from 5-decays into 
charmed particles. What can be measured is an integral over the endpoint 
region, 

dT(B^X u £u) 



dE f 



E 



dE P 



(54) 



where Eq is above the kinematic endpoint for B — > D Iv transitions, i.e. 
Eo > 2.3 GeV. The traditional way to extract a value of V u b from such 
a measurement is to compare T U (E ) with the predictions of various quark 
models. Since the endpoint region of the lepton spectrum is strongly affected 
by nonperturbative effects, this procedure suffers from a considerable amount 
of model-dependence. Currently, the value of V u b obtained following this 
strategy has a theoretical uncertainty of at least a factor 2 |42, 43fl . 

Based on the results obtained in this paper and in Ref. [P^| , we propose 
a new strategy to extract V u b with little model-dependence. The idea is to 
use the second equation in fl53|) to relate the integral T U (E ) to a weighted 
integral over the photon spectrum in rare decays. Defining 



T,(E l ,) = ±[dE.,(E,-E,,)<l T{B ^ X ^ 



m B 



Eo 



(55) 



and using 


that 


1 v tb v* 


1 rN ~ / 






2 




2 




v cb 




v tb v t * s 





V c b I , we find the remarkable relation 



(56) 



where ?]qcd contains radiative corrections, which have so far been neglected in 
this paper. The Wilson coefficient Cj(m b ) can be calculated in perturbation 
theory. In the standard model, it is known to next-to-leading logarithmic 
T6| . Hence, from a measurement of the integrated quantities 



accuracy 
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T u (Eq) and T S (E ) one obtains a direct determination of |V^&/Vd>|. The 
fact that the right-hand side in fl56[) must be independent of E Q provides a 
constraint, which can help in the analysis of the data. By taking the ratio of 
the integrated decay rates, we are able to reduce hadronic uncertainties to 
the level of power corrections, the leading ones being of order A/m^ ~ 0.1, 
where we take A ~ 500 MeV as a typical low-energy scale of the strong 
interactions. It is thus not inconceivable that the theoretical uncertainties 
in fl56|) can be controlled to a level of, say, 10-30%, corresponding to an 
uncertainty in of 5-15%. This would be a major improvement over 

the present situation. 

To achieve such an accuracy, it is necessary to study in detail the QCD 
correction factor ^qcd in which arises when radiative corrections are 
included in the operator product expansion of the inclusive decay rates. We 
shall briefly discuss the qualitative structure of these corrections; a more 
complete treatment will be presented elsewhere In the free quark decay 
model, the one-loop radiative corrections to B — > X s 7 and B — > X u i v 

For the 



the one-loop radiative corrections to B 
decays have been investigated by several authors |T7|, £|5 



integrated quantities Ti(E ), they have the general structure (i = u or s) 



rr ton (E )oc(l-y )e(l-y ){l 



2a, 
3tt 



ln 2 (l - y ) + a ln(l - y ) + 6, 



O 



:i-i/o) ; 



(57) 



where y = 2E /rrib- The Sudakov-type double logarithms are universal and 
enter both quantities with the same coefficient. This statement is true to all 
orders in perturbation theory fi9|| . Note that the result for T s (y ) is still 
of this form when one takes into account the finite strange quark mass. This 
effect is known to modify the endpoint behavior of the spectrum in the very 
endpoint region 1 — yo ~ p = 0(l/ml) [0. However, since the integral in 



(1551) extends over a larger region of order 1/rrib, the corrections induced by 
p 7^ are subleading. They are of the same magnitude as terms of order 
(1 — yo) 2 , which we neglect in (|57D. 

Since integrations over y and k + commute, we can incorporate bound- 
state corrections by convoluting the parton model result (|57D with the struc- 
ture function f(k + ), as we did for the tree-level expressions in ( PB[ ) and 
Tf). When we then take the ratio of T U (E Q ) and T S (E ), the large double- 
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logarithmic corrections cancel. We find 



VQCD 



Here 



2a, 
3tt 



dk + f(k A 



[a u - a s ) lnr + (b u - b s 



(58) 



lnr 



tn b {y -l) 



i - y + — in 



a. ^+ 
2/o H 



(59) 



— In r > — In 1 



(!*+/(*;+) l-y„+ — 
m b (ya-l) 

is a nonperturbative parameter of order Aqcd/^6- Using the fact that the 
logarithm is a monotonic function over the range of integration, and that 
1 — 2/o + A/mb = l— 2E /rriB up to corrections of order l/m|, we find the 
bound 

^) • (60) 

A precise calculation of lnr, however, requires some knowledge of the non- 
perturbative structure function f(k+). The coefficients a s and b s in (|58|) can 
be extracted from the analysis of Ali and Greub [[17]]. We obtain a s = 3/2 and 
b s = 2n 2 /3 + 1. In b s , we have neglected small contributions from operators 
other than O-j in the short-distance expansion of the effective Hamiltonian 
([I]). Unfortunately, the existing analytical calculations of the radiative cor- 
rections to the lepton spectrum in B — > X u £ v decays disagree on the value 



of b u . We find a u = 19/6 and b u 



7T 



23/12 + 6, where S 



to Jezabek and Kiihn while 5 ~ 3/8 according to Corbo [46 
correction factor //qcd, we find 



according 
For the 



^QCD 



2a, 
9tt 



35\ 



blar + ir + S - 



1 + 



1.54 



1.11 In — 
,0.2 



4 J 
0.22 5 



a, 
7T 



(61) 



Since the controversial quantity 5 enters with a small coefficient, the corre- 
sponding uncertainty in ^qcd is small. Using a s /ir ~ 0.08 and r ~ 0.2, we 
expect tjqcd — 1.12. We conclude that the perturbative corrections are not 
unexpectedly large, but they deserve further investigation. In particular, the 
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nature of the single-logarithmic corrections should be clarified. It is tempt- 
ing to interpret the lnr-term in ( |58D as a renormalization-group logarithm 
arising from the scaling from large scales /i 2 ~ m 2 down to smaller scales 
jj 2 ~ m 2 (l — y ) ~ m b AQCD, which are characteristic of the endpoint re- 
gion. If this interpretation is correct, it should be possible to resum these 
logarithms using renormalization-group techniques. At present, however, we 
cannot prove this assertion. 

Let us finally compare our new strategy to alternative model-independent 
determinations of V u b from exclusive S-decays. Using heavy quark flavor 
symmetry, one can in principle extract V u b from a comparison of exclusive 
decays B — > hiv and D — > h£u, where h is a light hadron such as tt 
or p (see, e.g., Refs. EO, BT1 152]). The problem is that the comparison 



must be done close to the zero recoil limit. Thus, one is restricted to a 
small fraction of phase-space. Presently, there is no convincing experimental 
evidence for exclusive charmless -B-decays. But even if such an analysis 
becomes feasible as high-statistics data from a 5-factory become available, 
the theoretical uncertainties will be of order l/m c instead of 1/mb- We thus 
believe that our new approach, based on a comparison of inclusive decay 
spectra, is preferable from the point of view of both theoretical uncertainty 
and experimental feasibility. 



7 Summary and Conclusions 

We have presented a systematic, QCD-based analysis of bound-state cor- 
rections to the photon spectrum in inclusive B — > X s 7 decays. Using the 
operator product expansion and the heavy quark effective theory, we are 
able to resum the leading nonperturbative effects into a universal structure 
function f(k + ), which describes the distribution of the light-cone residual 
momentum of the 6-quark inside the 5-meson. This formalism provides the 
generalization of the phenomenological concept of the "Fermi motion" in the 
context of QCD. We find that the moments of the structure function are given 
by a set of universal forward matrix elements of higher-dimension operators. 
The characteristic width of the photon spectrum is related to the expecta- 
tion value of the kinetic energy of the heavy quark inside the 5-meson. As 
a by-product, we obtain the bound Ai < for one of the parameters of the 
heavy quark effective theory. 
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The formalism presented here is rather general. So far, it has been ap- 
plied to the photon spectrum in rare decays and to the lepton spectrum 
in B — > X u lv transitions P5| . Applications to other processes such as 
B — > X c iv |53| or purely hadronic decays are possible, too. The fact that 



the leading nonperturbative corrections to inclusive decay spectra can be 
traced back to a universal structure function leads to interesting relations 
between different processes. In Sect. || we have shown that a weighted inte- 
gral over the photon spectrum in B — > X s 7 decays is related to an integral 
over the endpoint region of the lepton spectrum in B — > X u £ v decays. Based 
on this connection, we have proposed a model-independent way to extract 
the ratio |V^b/Vd,| of elements of the Kobayashi-Maskawa matrix. Hadronic 
uncertainties enter this determination only at the level of l/ra& corrections. 
We estimate that using this method one could extract V u b with a theoretical 
uncertainty of about 20%, which is an order of magnitude better than the 
present theoretical uncertainty in this parameter. 

Since the 6-quark structure function is a rather fundamental object, one 
should try to calculate it using nonperturbative techniques such as lattice 
gauge theory or QCD sum rules. We believe that the relations derived in 
Sect. [| will be helpful for such calculations. In particular, it may be of 
advantage to calculate the Fourier transform of the structure function, which 
is given by the forward matrix element of a bilocal operator. Any theoretical 
insight into the behavior of the structure function and its moments will have a 
direct impact on the analysis of inclusive decay spectra, from which one hopes 
to extract accurate values for some standard model parameters or, in the case 
of rare decays, even indications for new physics beyond the standard model. 
In this context, we emphasize that the phenomenological model of Altarelli 
et al. [pjj, which has been widely used to account for bound-state corrections, 
corresponds to a simple Gaussian model for the structure function and is not 
entirely consistent with QCD. A measurement of some of the moments of 
the shape function in semileptonic decays |2Sj , or of the photon spectrum in 
rare decays, could help to find deviations from this model. More generally, 
such a measurement would provide us with some fundamental QCD matrix 
elements and is thus interesting in its own right. 

Finally, we mention that before the formalism developed here can be 
applied to the analysis of data, it is necessary to include QCD radiative cor- 
rections. This can be done by calculating virtual and real gluon corrections 
in the free quark decay model, and convoluting the result with the struc- 
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ture function f(k+). This procedure has its subtleties, however, due to the 
presence of endpoint singularities in the perturbative expansion. We have 
indicated this for the important case of the ratio of integrated decay rates in 
(|56|) . A more complete discussion will be given elsewhere p4| . 
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